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Abstract. We explain how scale dependent renormalized quantities can be com- 
puted using lattice QCD. Two examples are used: the running coupling and quark 

J"""""- ' masses. A reliable computation of the yl-parameter in the quenched approximation 

ON ■ is presented. 

Ov 
q \ 1 Introduction 

[\ \ In a perturbative treatment, one takes as parameters of QCD the running 

coupling and running quark masses in a specific renormalization scheme. If 

we adopt a mass-independent scheme, the coupling and quark masses satisfy 

^—1 ■ rather simple renormalization group equations. Their asymptotic high energy 

^ \ behavior is given by the yl-parameter and the renormalization group invariant 

masses {Mf} (/ labeling the flavors). These parameters have the advantage 
that Mf are independent of the renormalization scheme and the /1-parameters 
of different schemes can be related exactly by 1-loop calculations. 

Once one considers QCD on the non-perturbative level, it is natural to 

f**«« ■ take N{ + 1 hadronic observables to be the parameters that fix the theory 

ON ' instead of A and {Mf}. We call this a hadronic scheme (HS). For QCD with u 

ff ■ and d taken as mass-degenerate and an additional s-quark, one may take for 

<**> ij example the proton mass, m p , the pion mass, m 7!l and the Kaon mass, mK, 

Jv, ■ as the basic parameters. Of course, the two sets of parameters are related 

g_) ' once a non-perturbative solution of QCD is available. It is a challenge for 

^^ lattice QCD to provide this relation with completely controlled errors. In 

^. ' this brief report we describe a method that allows to solve this problem. The 

• *h . numerical results (sect. 4) are for quenched QCD, i.e. the fermion determinant 

/\ ' is omitted in the path integral; one may think of this as a version of QCD 

3 • with JVf = flavors of sea-quarks. 

The problem described above is the connection of the short distance and 
long distance dynamics of QCD. As such it appears at first sight to be in- 
tractable by numerical simulations, since it requires that various scales are 
treated on one and the same lattice. This is impossible, if discretization 
errors are supposed to be under control! However, one may define a running 
coupling and running quark masses in QCD in a finite volume of linear size L 
with suitably chosen boundary conditions (cf. sect. 3). Then the masses and 
coupling run with an energy scale q — l/L and their evolution can be com- 
puted recursively. The only requirement to keep discretization errors small in 
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each step of the recursion is L ;§> a which is easy to satisfy (a denotes the 
lattice spacing) |y, || . 



2 Strategy 

We start by giving an overview of the strategy to compute short distance 
parameters in figure 1. One first renormalizes QCD replacing the bare pa- 

L max = O(ifm) : HS — » SF(g = 1/L max ) 

i 
SF(q = 2/£ max ) 

I 



i 

SF(q = 2™/L max ) 
PT: | 

PT 

jet - physics < — vIqcd , M 

Fig. 1. The strategy for a non-perturbative computation of short distance param- 
eters. 



rameters by hadronic observables. This defines the HS introduced above. It 
can be related to the finite volume scheme (denoted by SF) at a low energy 
scale q = l/£ max , where L max is of the order of I fm. Within this scheme one 
then computes the scale evolution up to a desired energy q = 2"/L max . It is 
no problem to choose the number of steps n large enough to be sure that one 
has entered the perturbative regime. There perturbation theory (PT) is used 
to evolve further to infinite energy and compute A and {Mf}. Inserted into 
perturbative expressions, these parameters provide predictions for jet cross 
sections or other high energy observables. In figure f , all arrows correspond to 
relations in continuum QCD; the whole strategy is designed such that lattice 
calculations for these relations can be extrapolated to the continuum limit. 

For the practical success of the approach, the finite volume coupling and 
quark masses must satisfy a number of criteria. 

— They should have an easy perturbative expansion, such that the /3-function 
(and r-function, which describes the evolution of the running masses) can 
be computed to sufficient order in the coupling. 

— They should be easy to calculate in MC (small variance!). 

— Lattice effects must be small to allow for safe extrapolations a/L — > 0. 
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Careful consideration of the above points led to the introduction of the run- 
ning coupling, g, and quark mass, to, through the Schrodingcr functional 
(SF) of QCD |,|,|,1. 

3 Schrodinger functional scheme 

We describe the definition of g and to in a formal continuum formulation. The 
SF is the partition function of the Euclidean path integral on a hyper-cylinder 
with Dirichlet boundary conditions in time, (P± = ^(1 ± 70)) 

A k (x)=Ck(x), P+i/)(x) = p(x), jp_(x)P- = p(x) atx =0, ,, 
A k {x) = C£(x), P^{x) = //(x), i>(x)P+ = p'(x) at x = L . [ ' 

Here, Af., . . . , p' are classical prescribed boundary fields. In space, the gauge 
fields are taken periodic under shifts x — > x + Lk and the fermion fields are 
periodic up to a phase 9. 

The renormalized coupling is defined through the response of the SF to an 
infinitesimal change of the boundary gauge fields Ck,C' k . Since the boundary 
fields are taken with a strength proportional to 1/L, there is no scale apart 
from L and the coupling g(L) runs with L. Details have been discussed in 
several publications ||, ||] . 

A natural starting point for the definition of a renormalized quark mass 
is the PCAC relation, 

d tl Al{x) = 2mP a (x) , (2) 

which expresses the proportionality of the divergence of the axial current, 
A 1 1, to the pseudo-scalar density P a , where 

P a (x) =^(x) l5 ±r a iP(x) 7 Al=^(xh^T a ^(x) , (3) 

and r a are the Pauli-matrices acting in flavor space (we take N{ = 2 degen- 
erate flavors from now on). Renormalizing the operators in equation |2| 

(A R ); = Z A A* , P£ = Z P P a , (4) 

we define a renormalized quark mass by 

to = ——to . (5) 

z P 

Here, to is to be taken from equation inserted into an arbitrary correlation 
function. It does not depend on the chosen correlation function, since equation 
|2| is an operator identity. The scale independent renormalization Za can be 
fixed through current algebra relations (also in the lattice regularization || [7| ) 
and we are left to give a normalization condition for the pseudo-scalar density, 
Pr,. The running mass, to, then inherits its scheme- and scale-dependence 
from the corresponding dependence of Pr. 



R. Sommer 





space 



Fig. 2. /p (left) and /i (right) in terms of quark propagators. 



We define Zp in terms of correlation functions in the SF. To start, we 
construct (isovector) pseudo-scalar fields at the boundary of the SF a, 



O a 



d 3 i 



d 3 v C(u) 75 ir a C(v) 



(6) 



C, C Their precise definition S involves a 
. After taking this functional derivative 



from the "boundary quark fields 

functional derivative, e.g. ((x) = j-— T 

one sets the boundary values, p... ,to zero. Analogously one defines a field 0' a 

which resides at xq = L. These fields are used in the correlation functions 



fp(x ) = -UP a (x)O a ) , h = (0' a O a ) , 



(see figure 2). In the ratio 



Const . V7l //P (x) \x =L/2 



(7) 



(8) 



Fur- 
The 
This 



the renormalization of the boundary quark fields H cancels out; it can there- 
fore be taken as the definition of the renormalization constant of P. The 
proportionality constant is to be chosen such that Zp = 1 at tree level, 
ther details of the definition are C = C = and a specific choice for 9 
renormalization constant is to be evaluated for zero quark mass, m. 
defines a mass- independent scheme. 

By construction, the SF scheme is non-perturbative and independent of a 
specific regularization. For a concrete non-perturbative computation, we do, 
however, need to evaluate the expectation values by a MC-simulation of the 
corresponding lattice theory. For the details of the lattice formulation we refer 
to B but mention that it is essential to use an 0(a)-improved formulation in 
order to keep lattice spacing effects small M . We now explain the computation 
of the scale dependence, omitting the matching to the HS scheme [pi pf for 
lack of space. 



At this point, the SF provides us with an important advantage compared to other 
boundary conditions, namely we can project the boundary quark fields onto zero 
momentum. As a result, the correlation functions vary slowly with xo, leading to 
both small statistical and small discretization errors. 
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4 Scale evolution, A and RGI masses 

Each step in the recursive computation of the scale evolution consists of: 

1. Choose a lattice with L/a points in each direction. 

2. Tune the bare coupling, go, such that the renormalizcd coupling g 2 (L) 
has the value u and tune the bare mass, mo, such that the mass (equation 
g) vanishes; compute Zp. 

3. At the same values of go, mo, simulate a lattice with twice the linear size; 
compute u' = g 2 (2L) and Z P . This determines the lattice step scaling 
functions for the coupling, S(u, a/L)=u', and for the mass, Sp(u, a/L) = 
Zp/Zp. 

4. Repeat steps 1.-3. with different resolutions L/a and extrapolate cr(u) = 
Unia/L^o ^{ u i O'/L), <Jp{u) = lim a /£^o Sp{u, a/L). This extrapolation 
can be done by allowing for small linear a/L-discrctization errors H . 
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Fig. 3. Non-perturbative running coupling and quark mass as a function of q = 1/L. 
Uncertainties are smaller than the size of the symbols. 



In this way, u(u), 1.1 < u < 2.4 and ap(u), 1.1 < u < 3.48 have been 
computed in the quenched approximation. The scale evolution is monotonous 
in this range: <j(u) > u, op{u) > 1. a and up can therefore be used to follow 
figure 1 downwards (up in the energy scale). 

Starting from the box size i m ax, defined by g 2 (L max ) = 3.48, we use 



g 2 {Lk-i)=v{g 2 (L k )), m(Lk-i) = (Tp(g 2 (L k ))m(L k ), L k 
to compute the sequences j|) 

g 2 {L k ), m(Lfe)/m(i max ), k = 0, 1, ...8. 



2~ k L r 



(9) 



(10) 



In figure 3, these are compared to the perturbative evolutions (starting at 
the weakest coupling) . It is surprising that the perturbative evolution is quite 



6 R. Sommer 

precise down to rather low energy scales. This property may of course not be 
generalized to other schemes but has to be checked in each scheme separately. 
On the other hand, the figure clearly demonstrates that the SF-scheme is 
perturbative in the large energy regime, say q > 50GeV. We may therefore 
start from the exact relations 

A = q (b - 9 r bl/ ^e-W>^exp\- f 9 dx\d* + ^--k\) <") 



M = fn (2b g 2 )- do/2b ° exp j - f dg[j 



o 
(a) do_] 

Ws) b g\ 



fl(x) bgx 3 bf^c 



(12) 



insert for the /3-function (r-function) its 3-loop (2-loop) approximation and 
set q to the maximal value reached, to compute A and fn/M. One can eas- 
ily convince oneself that the left over perturbative uncertainty is negligible 
compared to the accumulated statistical errors. 

After converting to the MS-scheme one arrives at the result M 

44 = 251 ± 21MeV , (13) 

MS ' v ' 

where the label '°' reminds us that this number was obtained with zero quark 
flavors. The error in equation [L3| accounts for all uncertainties including the 
extrapolations to the continuum limit that were done in the intermediate 
steps. 

This talk is based on research done in the framework of the ALPHA 
collaboration |3]. I thank in particular S. Capitani, M. Guagnelli, M. Luscher, 
S. Sint, P. Weisz and H. Wittig for the joint work. 
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